Charge and Spin Transport at the Quantum Hall Edge of Graphene 
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Landau level bending near the edge of graphene, described using 2d Dirac equation, provides 
a microscopic framework for understanding the quantum Hall Effect (QHE) in this material. We 
review properties of the QHE edge states in graphene, with emphasis on the novel phenomena 
that arise due to Dirac character of electronic states. A method of mapping out the dispersion of 
edge states using scanning tunneling probes is proposed. The Zeeman splitting of Landau levels is 
shown to create a particularly interesting situation around the Dirac point, where it gives rise to 
counter-circulating modes with opposite spin. These chiral spin modes lead to a rich variety of spin 
transport phenomena, including spin Hall effect, spin filtering and injection, and electric detection 
of spin current. The estimated Zeeman spin gap, enhanced by exchange, of a few hundred Kelvin, 
makes graphene an attractive system for spintronics. Comparison to recent transport measurements 
near v — is presented. 
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INTRODUCTION 

Isolation and gating [l[ of graphene, a monolayer of 
graphite, has enabled the observation of interesting trans- 
port effects, resulting from Dirac fcrmion-likc character 
of excitations. In particular, graphene hosts an integer 
quantum Hall effect (QHE) with unusual plateau struc- 
ture [2j, [3j • It was found that the QHE plateaus in mono- 
layer are arranged symmetrically around the neutrality 
point, occurring at filling factors which are half-integer 
multiples of four, which is the combined spin and valley 
degeneracy of graphite. This form of QHE came to be 
known as the anomalous or half-integer QHE. 

The simplest framework allowing to understand this 
behavior of QHE is provided by the structure of Landau 
levels (LL) of a 2d massless Dirac equation [28] which has 
the particle-hole symmetric spectrum 



eo = Hvoy/2eB/hc 



(1) 



with the sign plus for positive n, and minus for negative 
n. Here B is the magnetic field and vq w 8 x I0 7 cm/s 
is the velocity at the graphite Dirac point K or K' . 
Due to spin and valley degeneracy (we shall discuss 
the role of Zeeman splitting below), each of the levels 
(JTJ contributes a step of 4e 2 //i to the quantized Hall 
conductivity. The particle-hole symmetry e — ► — e of 
the spectrum (JTJ), with the n — level positioned at 
e = 0, suggests that the QHE plateaus must occur at 
v = ±2(2n + 1) = ±2, ±6, ±10, which is indeed what 
is observed in experiment 0, Q • The anomalous QHE in 
graphene can be understood in a more fundamental way 
in terms of a quantum anomaly of the zeroth Landau 



[28] The result JTJ follows from the well-known representation of the 
Schrodinger-Pauli Hamiltonian as a square of a massless Dirac 

Hamiltonian, -HgP = ("^(Pi ~ t^*)) 1 ^ ne Dirac spectrum 
(TTJ is then obtained by taking a square root of the nonrelativistic 
Landau level spectrum. 



level [4| . The special character of the monolayer spec- 
trum ([1]) is underscored by the difference between the 
QHE properties observed in graphene monolayer and bi- 
layer systems [5[ . 

One of the most dramatic consequences of the Dirac 
LL spectrum (JTJ is the appearance of the new energy 
scale £o- The square root dependence £o oc \[B leads 
to a much larger level spacing than that for electrons 
with the quadratic dispersion conventional for semicon- 
ductors. For typical magnetic field of 10 T, the separa- 
tion of the lowest LL (n = 0, 1) is quite large, AE = 
£o ~ 1000 K, which enables QHE to persist up to room 
temperature [6( . 

To gain insight into the microscopic origin of the 
anomalous QHE, it is useful to develop the edge-states 
approach, which provides an intuitive and simple pic- 
ture of the conventional QHE [3] ■ The edge states for 
graphene were studied using a numerical solution of the 
tight-binding model § and also with the help of the Dirac 
equation (il. floj| . It was found that the energy levels (JTJ) , 
valley-degenerate in the graphene bulk, are split near the 
edge due to valley mixing at the boundary. Interestingly, 
it turns out that the structure and dispersion of the edge 
states depend on particular crystallographic orientation 
of the edge. For the so-called armchair edge, a simple 
particle- hole symmetric splitting was found [9( . Particle- 
like and hole-like states have different sign of energy dis- 
persion, giving rise to counter-propagating modes with 
opposite chirality. A somewhat more complicated situa- 
tion occurs near the zigzag edge, where dispersing edge 
states coexist with an additional dispersionless surface 
state [§]. Despite these differences, the armchair, zigzag 
and other edges have same numbers of dispersing edge 
states of both chiralities. As we discuss below in Sec. II, 
this ensures the universal half-integer character of QHE 
in graphene. 

One unique aspect of QHE in graphene is that its elec- 
tronic states, owing to the monolayer character of this 
material, are fully exposed and, similar to the surface 
states of 3d materials, can be investigated by scanning 
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tunneling microscopy (STM) probes with atomic reso- 
lution. Moreover, some useful information can even be 
obtained by imaging the top layer of 3d graphite, as 
demonstrated by recent STM studies of Dirac Landau 
levels 



edges 




and of electron states near atomically thin 
141 ] . In Sec. II we discuss new possibilities for 
scanning experiments that arise in graphene. The charac- 
teristic spatial scale of the states in a Landau level, given 

1 /2 

by the magnetic length £b = (hc/eB) , is about 8nm 
for the field of 10 T. Being large compared to the STM 
spatial resolution, it allows to image individual electronic 
states with sub-is resolution and, in particular, to study 
Landau level bending near graphene edge. This bending 
in fact mimics the edge states momentum dispersion, due 
to the position-momentum duality of the Landau levels. 
As we shall see, the STM technique has sufficient resolu- 
tion to map out the dispersion of QHE edge states. 

Another novel feature of graphene is the simultaneous 
presence of the QHE edge modes of opposite chiralities, 
propagating in the opposite directions. Being particle- 
like and hole-like, they occur at the energies e > and 
e < 0, respectively. For electron density detuned from 
the neutrality point, v = 0, only one of the chiralities 
contributes to transport. However, as we shall see, near 
v = the states of both chiralities can participate in 
transport, leading to rather unusual transport properties. 
In particular, in the presence of Zeeman spin splitting of 
Landau levels, the v = state features an energy gap 
in the bulk and, simultaneously, a pair of edge states of 
opposite chirality and opposite spin polarization [^, UH . 
These states carry spin-up and spin-down electrons in 
the opposite directions along graphene edge, exhibiting 
quantized spin Hall effect but no charge Hall effect, owing 
to the particle-hole symmetry at v = 0. As we discuss in 
Sec. Ill, these spin-polarized chiral edge states exhibit in- 
teresting spin transport phenomena, such as spin filtering 
and spin injection, whereas the spin Hall effect provides 
a natural tool for the detection of spin current. 

Interestingly, the counter-propagating edge states 
manifest themselves directly in charge transport. As we 
discuss in Sec. IV, near v — longitudinal resistivity re- 
mains finite, p xx > h/e 2 , while the Hall effect, which 
is absent at the particle- hole symmetry point v = 0, 
appears at nonzero v due to conductivity in the bulk. 
This leads to Hall resistance p xy changing sign at v = 
without exhibiting a clear plateau. The bulk conduc- 
tivity short-circuits the edge transport and suppresses 
longitudinal resistivity, leading to a prominent peak in 
p xx near v = and a plateau in a xy These predictions, 
as well as the behavior of resistance fluctuations, which 
are enhanced near zero v, are in agreement with experi- 
ment [3, [nl • 

Looking beyond graphene, interesting massless Dirac 
fermion states have been predicted a while ago [HI at in- 
terfaces of narrow-gapped HgTe and PbTe semiconduc- 
tors. In Ref. pjl], which anticipated many of the features 



of electronic states in graphene, 2d Dirac states occur in 
a band-inverted heterojunction plane due to spin orbit 
interaction in 3d bulk, in the absence of magnetic field. 

We note also that there have been interesting predic- 
tions of quantized spin Hall effect in certain 2d insula- 
tors In particular, the proposal of Ref.[20j can 
be viewed as a 2d version of the 3d situation discussed 
in Ref. 18]. In these schemes an energy gap forms in 



the bulk due to spin-orbit even in the absence of mag- 
netic field, similar to Ref. 18]. At the same time, counter- 
propagating modes carrying opposite spins appear on the 
edge, which are responsible for the quantized spin Hall 
current. As was emphasized in Ref. 22], there are gen- 
eral symmetry requirements, rooted in the time-reversal 
symmetry, protecting counter-propagating gapless exci- 
tations at an insulator boundary. In particular, certain 
Zi invariants must exist, which are realized as spin <j z 
component in the case of graphene [UQjl], and were linked 
to supersymmetry for heterojunction systems [lit ]. 

We note that the nature of the edge states of Refs. 
liol ] is essentially identical to ours. The only differ- 
ence, apart from different size of the energy gap, is that 
Rashba spin-orbit may cause backscattering in our case 
(see Sec. Ill) but not in the situation of Refs. ID, ljj 2p| . 
while magnetic impurities will cause backscattering in 
both cases. In particular, our discussion in Sec. IV of p xx 
and p xy and our proposals to operate and detect spin cur- 
rent can serve as diagnostic tools should these schemes 
become realized experimentally. 

The rest of the paper is organized as follows. In Sec. II 
we introduce the edge states using the Dirac equation 
framework. We focus on the two main edge types, arm- 
chair and zigzag, however we emphasize the generic as- 
pects that are applicable to other edges. Then we discuss 
the anomalous QH effect and the possibility of imaging 
the edge states with scanning tunneling probes. In Sec. Ill 
we focus on spin-polarized chiral edge states and related 
spin transport phenomena. We also comment on spin re- 
laxation mechanisms and present estimates of the spin re- 
laxation time. In Sec. IV we introduce a transport model 
which accounts for both edge and bulk transport. This 
allows us to connect the chiral spin-polarized edge picture 
with recent transport measurements near the neutrality 
point [lEll3- 



II. DIRAC QHE EDGE STATES 

Here we analyze electron states near zigzag and arm- 
chair edges (Figfjji), the two most common graphene 



edge types, using massless Dirac model [21|. This exer- 



cise, which amounts to setting the boundary conditions 
for the Dirac spinor and solving an appropriate Id eigen- 
value problem, provides a fully microscopic picture of 
graphene QHE. Below we use this approach to illustrate 
the interplay between the QHE edge states and surface 



3 



(a) Armchair edge (b) 




FIG. 1: (a) Graphene lattice with armchair and zigzag edges. 
Sublattices A and B are marked, (b) Graphene hexagonal 
Brillouin zone with Dirac valleys K and K' . (c) The two 
linearly independent zero-energy Bloch functions for the K 
valley used in J2j, where r = e 2?rl/3 . For the K' point, the 
zero-energy Bloch functions are obtained from those shown 
by complex conjugation. 




FIG. 2: (a) Electron energy spectrum near the armchair 
boundary obtained from the Dirac model, Eq.([2]). The bound- 
ary condition, Eq.©, lifts the K, K' degeneracy, forming 
particle-like and hole-like counter-propagating edge modes. 
The odd integer numbers of the modes yields the "half- 
integer" QHE. (b) Spin-split graphene edge states, propagat- 
ing in opposite directions at zero energy: the blue (red) curves 
represent the spin up (spin down) states. 



states for the zigzag edge, and discuss the possibility of 
imaging the edge states with scanning tunneling probes. 

Let us first recall how low-energy graphene excita- 
tions are obtained in the tight-binding model 2l| near 
the Dirac valleys K and K', located at the two non- 
equivalent Brillouin zone corners (see FigOJb)). There 
are two linearly independent zero-energy Bloch functions 
for each of the points K, K', each residing only on one 
sublattice (A or B) and vanishing on the other sublat- 
tice. Our choice of Bloch functions for the K valley is 
shown in FigOJc) ■ The Bloch functions for the valleys K 
and K' are related by complex conjugation. 

The wave function of low-lying excitations near K and 
K', is written as a superposition of these four zero-energy 
Bloch functions multiplied by slowly varying envelope 
functions uk, Vk, —uk>, —vk>, with u and v being 
the wave function amplitudes on the sublattice A and 
B. (Our choice of the signs for the K' valley is conve- 
nient for treating an armchair boundary, as we shall see 
below.) The envelope functions uk, vk, and Uk>, Vk' 
describe excitations near K and K', respectively. The 
effective low-energy Hamiltonian, obtained by keeping 
only lowest-order gradients of u and v, takes the massless 
Dirac form [29]: 



H 



K 



o p+ 
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H 



K> 



IV 





-P-\ 



V- 




, (2) 



where p±=p x ± ip y , p^^p^- {e/cjA^. 

Landau levels in an external -B-field, described with 
the gauge A x = —By, A y = 0, can be obtained for the 



states with the x dependence e WmX from Id Hamiltonians 



H K . 



K' 



v/2 







(y 
o 



y*) 



(3) 



1/2 

where So = Tlvq {2eB /tie) and j/„ = — p x . Here y and 
p x are measured in the units of Ib and %/Ib, respectively. 
The spectrum of Hk,k', Eq.Q, yields the Dirac Landau 
levels, Eq.([T]), where the eigenstates for the two valleys 
K(K') are given by 

(uK,n>VK,n) = A (c„<£„_i (y - ip n (y - y*)) , (4) 

{uK>,n,VK',n) = A (<p n (y - V*), Cntyn-X (y ~ V*)) -(5) 

Here ip n (z) are the eigenfunctions of the magnetic oscil- 
lator, n = 0, 1, the normalization factor A equals I for 
n = and l/v2 for n ^ 0, with Co = and c n ^o = L 
Note that zeroth LL states reside solely on B sublattice 
for valley K and on A sublattice for valley K' [30]. 

We now analyze how LL spectrum is modified near 
the armchair edge. We consider graphene sheet in the 
halfplane y < with an armchair edge parallel to the 
x axis (see FigOJa)). Energy levels near the edge are 
determined from the Dirac eigenvalue equations Eip = 
HjcK'ip-i where ip — (u,v), and Hk,k' are given by 
Eq.([3]). To analyze this eigenvalue problem, we exclude v 
components and consider eigenvalue equations with spec- 
tral parameter A = (E/eq) 2 for u components: 



1) U K = Xu K , 



2 v -9y + (y - y* 
| + {y- y«) 2 - 1) u K , = \u K , 



(6) 



The boundary conditions for Eqs.© can be obtained 
from the tight-binding model, which is valid up to the 



[29] We note that the Hamiltonian used in our paper [j| is related to [30] This property is specific for the zeroth LL, making the splitting 
the Hamiltonian JJJ by variable transformation uj^ jf' U K,K' > °f the n = LL due to Coulomb interaction distinctly different 

V K.K' ~* ~i v K K' ■ from that of other LLs (see [27l. |28H ). 
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very last row near y = 0, by setting the wave function 
equal zero at the boundary. Since the armchair edge has 
lattice sites of both A and B type (see FigfjTa)), the wave 
function on both sublattices should vanish at the edge. 
In terms of the envelope functions uk,k' , uk,K' , taken at 
y = 0, this condition translates into 



UK = U K > 



VK = V K > ■ 



(7) 



We obtain a pair of differential equations ^ on the semi- 
axis y < 0, coupled at the boundary via Eq.([7]). To 
simplify this problem, let us map Eq.© for uk 1 onto the 
positive semi-axis y > 0, by UK'{y) — > UK'{—y), while 
keeping uk on the negative semi-axis y < 0, and treat it 
as an eigenvalue problem in the domain — oo < y < oo 
with the wavefunction given by uk at negative y and by 
UK' at positive y. The first boundary condition, Eq.([7]), 
then means that the wavefunction is continuous at y = 
0, while the second condition implies continuity of the 
derivative du/dy. (This can be seen by expressing Vjc,K' 
in terms of uk,k> using Eqs.Q.) Thus we obtain a Id 
Schrodinger problem in the potential [j| 



V{y) = ^{\y\+y*) 2 -~sgn{y), 



(8) 



defined on the entire y axis. After finding the spec- 
trum A(y*) numerically, we obtain the energy levels of 
the Dirac fermions as 



E(p x ) = ±£a\ZMjJ*j, y* = -p x , 



(9) 



whereby the particle-hole symmetry is restored due to 
the two possible signs of E{p x ). The dispersion J9J is 
illustrated in Fig(2k- We note that the LL double valley 
degeneracy in the bulk is lifted at the boundary. 

The particle-hole symmetric edge states spectrum in 
Figl2k instantly explains the "half-integer" Hall quanti- 
zation in graphene. Indeed, for any electron density with 
integer v in the bulk there is an odd number of the edge 
modes crossing the Fermi level, which means that the 
Hall conductivity is quantized as a xy — 2{2n + l)e 2 /h, 
where the factor two accounts for spin degeracy. 

Spin degeracy of the Landau levels is lifted by the Zee- 
man interaction, which is substantial in graphene, 



E z =gp B B^50K, g « 2, 



(10) 



for B = 30 T (compare to g ~ 0.1 in GaAs quantum 
wells). Zeeman-split edge states, depicted in Fig|2b, have 
interesting characteristics for electron density near neu- 
trality, v = 0. At this density the state in the bulk is 
spin-polarized, with the Zeeman gap further enhanced 
by exchange (see below). The two branches of counter- 
propagating edge states near e = 0, carrying opposite 
spin, have interesting properties that will be discussed in 
more detail in Sec. III. 

We now analyze the zigzag edge, which even at B = 
hosts a band of dispersionless zero-energy states bound 



to the edge 23]. We shall refer to these states as surface 
states, to distinguish them from the dispersing QHE edge 
states. As we shall see, the surface states contribute to 
the splitting of n — LL near the zigzag edge, in agree- 
ment with the tight-binding calculations [|| . 

We consider graphene sheet in the half-plane x > 0, 
with its first row consisting of B atoms (see FigQJi). For 
the states with the y dependence e lPvV , with A x = 0, 
A v = Bx, from ([2|) we obtain Id Hamiltonians 



H 



K,K' 



SO 

V2 



d x ± (x — x* ) 

-d x ±{x — £*) 



(11) 



where x„ = p y . Similarly to the armchair case, the 
spectrum can be found from the eigenvalue equation 
Eip = Hk,k'^, where ip = (u,v), which should be sup- 
plemented with the boundary conditions. For our zigzag 
edge the wavefunction have must vanish on all A sites at 
x = 0. For that both envelope functions uk, uk 1 have to 
vanish at the boundary, 



11 k = 0, Uk> 



0. 



(12) 



Excluding v components, we obtain two separate eigen- 
value problems for the spectral parameter A = (E/eo) 2 , 



| (-d x + (x - x*) 2 + l) uk = Auk, 
3 {-dx + {x- a;*) 2 - l) u K > = \u K >, 



(13) 



where both uk and uk' satisfy the hard wall boundary 
conditions (fl~2|) . The amplitudes vk,k> on the B sub- 
lattice can be expressed via ampliudes uk.k 1 on the A 
sublattice and eigenenergy E, 



vk = (so/V2E) {-d x + (x- X*)) U K , 
uk' = (eo/V2E) {-dx - {x - x* )) u K >- 



(14) 



The eigenvalue problems (|13|) with the hard-wall bound- 
ary conditions (| 1 2[) are familiar from the theory of edge 
states in the conventional QHE Q , and their spectrum 
A(x*) can be found numerically. The Dirac fermion en- 
ergy dispersion E{p y ) = ±eoy/\{x*) is shown in Fig[3] 

The behavior of n ^ LL's is similar to the armchair 
case: there are two branches of the edge states, one for 
each valley, degenerate in the bulk, x* 3> 1, which split 
near the edge. The zeroth LL, however, coexists with the 
surface state, which makes its behavior rather peculiar 
and different for the two valleys. 

In the K valley, which we discuss first, the zeroth LL K 
states reside solely on the B sublattice, see Eq.Q, and 
therefore automatically satisfy the boundary condition 
uk = 0. Thus there are zero-energy states for arbitrary 
values of x*, of the form ipa{x — a;,) cx e -( x ~ x *) I 2 . Let 
us consider the states with x„ far outside the graphene 
half- plane, x* <C — 1. Not too far from the boundary, 
such states can be approximated by an exponential 



v K {0 < x < \xj) oc e |a *l x , u K {x) = 0. 



(15) 
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Sublattice B Sublatticc A 




FIG. 3: Electron energy dispersion for the zigzag boundary 
obtained from the Dirac model, Eq.([2}, (|12|l , (a) Spectrum 
for the K valley. The zeroth LL morphs into dispersionless 
surface mode [23| near the edge, (b) Spectrum for the K 1 
valley. The zeroth LL mixes with the surface mode, giving 
rise to two branches of dispersing QHE edge states. 

which is identical to surface state wave function [23j ]. 
Thus the zeroth LL for valley K near the edge transforms 
into the surface mode. Being dispersionless, this mode 
does not contribute to the edge transport. The edge state 
spectrum for the valley K is displayed in FigGJi. 

Now let us consider the zeroth LL for the valley K'. 
For a;* 3> 1, we approximate the ground state of the 
oscillator (|13l) with the hard-wall boundary condition as 

uk'{x) = tpx, (x) « (p (x - x*) - tp (x + £*). (16) 

The ground state energy Ao(x*) is then approximated by 

Mx*) « (h), (17) 

where h = j(—d 2 + (x — x*) 2 — 1) is the effective Hamil- 
tonian for uk> component, Eq. (|13[) . and (...) denotes av- 
eraging over the normalized wave function (|16p . 

Evaluating (h) for x* ^> 1, when the state (|16[) has 
unit norm with exponential accuracy, we obtain 

A (x*) xsx*n- 1/2 e- x *. (18) 

From the relation E = ±£qvA, we find the energies for 
the two branches of dispersing edge states 

E ± (x.) « ±(2a;») 1 / 2 7r- 1 / 4 e^/ 2 £ . (19) 

Plugging this expression in Eq. (|14[) , we obtain the wave 
function on B sublattice for these two branches, 

v K , = ±xl /2 e~ x ' x . (20) 

which is again the surface state wave function (compare 
to Eq. (fT5)) ). We therefore conclude that for the K' valley 
the zeroth Landau level and the surface state mix giving 
rise to two dispersing edge modes. This is in agreement 
with the spectrum displayed in Fig|3jD. We see that, al- 
though the Dirac model is applicable only in a small part 
of the Brilloin zone, near points K and K' , it provides 
a description of the states at the zigzag edge, including 



FIG. 4: Position-dependent tunneling spectroscopy of 
graphene near the zigzag edge for sublattices A and B. Due 
to the momentum-position duality, the dependence of LDOS, 
Eqs. (|21[l . on the the distance to the edge x mimics the mo- 
mentum dispersion of the edge states. Note the difference of 
the spectra for the two sublattices near the edge and presence 
of surface state for sublattice B. 



the surface state, which is in agreement with the results 
of the tight-binding model of Ref. Q . The surface mode 
in the vicinity of K and K' is given by Eqs. (fT5|) ,(|2"0" ]) . 

Interestingly, the A and B sites contribute equally 
to the splitting of the zeroth LL, J x>0 \uk> \ 2 dx — 
f x>0 \vk' \ 2 dx. This is somewhat counterintuitive, since 
this LL is solely on the sublattice A in the bulk, while 
the surface mode is solely on the sublattice B. This 
equal participation property can be understood as fol- 
lows. The spinor states [uk 1 > i'RK"' ) with Uk>, vk' given 
by Eqs. (fl"6|) . (f20")) . are eigenstates of the Dirac Hamilto- 
nian with the boundary condition (fT2")l . with the energies 
±E. Thus these states are orthogonal, which implies that 
the integrals of |ux'| 2 and |fif| 2 are equal. We further 
note that the integral of the square of the B component 
of our edge state wave function ((20)) over x > indeed 
equals one, in agreement with our choice of normalization 
on the A sublattice in Eq. (fT6|) . 

To sum up, for the zigzag edge, the zeroth LL gives 
rise to two dispersing edge states for one of the valleys, 
while for the other valley the zeroth LL morphs into the 
dispersionless surface mode which does not contribute to 
the edge current. Therefore, despite the presence of the 
surface mode, the number of dispersing QHE edge states 
with e > and e < for the zigzag boundary is the same 
as for the armchair boundary, giving rise to "half-integer" 
quantization of Hall current. 

Finally, we briefly discuss how the edge states in 
graphene can be investigated using the STM tech- 
nique [H E3, El G3]. Due t0 the Landau level 
momentum-position duality relation, p y = (K/£g)x*, the 
edge state momentum dispersion shown in Figs trans- 
lates into the excitation energy dependence on the dis- 
tance from the edge. The characteristic scale for the 
latter is set by the magnetic length £g, which for typical 
fields is about 50-80 times greater than the spatial resolu- 
tion of STM instruments on graphite surface. This makes 
STM technique particularly convenient for this kind of 
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studies. 

A link between the edge states dispersion and the 
position-dependent tunneling spectroscopy can be estab- 
lished as follows. We shall use the solutions for the edge 
state wave function given above to calculate the local 
density of states (LDOS) near the zigzag edge (other 
edge types can be dealt with similarly). For each of the 
graphene sublattices LDOS is given by 



p A (E,x) =J2 a \u a (x)\ 2 6(E-E a ), 
p B (E,x) =E a K(x)\ 2 6(E - E a ), 



(21) 



where x is the distance from the edge, and a denotes 
the set of eigenstates of the K and K' Hamiltonians (fTTj) 
with the hard- wall boundary condition (fT2"j) . 

Using the eigenfunctions u a (x), v a (x) and the energies 
E a found from Eqs. (fl"3|) .(fl"4 |) as discussed above, we ob- 
tain LDOS for the A and B sublattices which is displayed 
in FiglU We see that the position-independent Landau 
level bands, dominating LDOS far from the edge, bend 
away from e — near the edge. This bending mimics the 
edge states momentum dispersion shown in Figl3l Note, 
however, that LDOS is nonzero only for x > 0, whereas 
the edge state momentum p y can be both positive and 
negative. The spatial width of the bending bands is de- 
termined by the width of the eigenfunctions u a (x), v a (x), 
which is of the magnetic length scale. 



III. SPIN-POLARIZED CHIRAL EDGE STATES 
AND SPIN TRANSPORT. 

As we noted above, at the neutrality point v = 
graphene hosts gapless spin-polarized edge states (see 
FigEtb)). The Zeeman energy gap in the bulk, Eq. (fT0|) . 
is enhanced by the Coulomb interaction. A Hartree-Fock 
estimate of this enhancement [9J gives 



A = 



(1 -a)ea « 0.456- (1 - a)e , 



(22) 



where k = 1 + ire 2 /2frvo ps 5.24 is RPA screening func- 
tion, and the parameter < a < 1 describes relative 
strength of Coulomb and exchange correlations. Assum- 
ing a = 0, i.e. ignoring correlations of electrons with op- 
posite spin, we obtain a spin gap A ~ 800 K for B = 30 T. 
Taking into account the substrate dielectric constant, 
e 2 — > j^j-e 2 , changes the result only slightly (^j- = 0.36 
for Si02). This approximation, while pointing at a cor- 
rect order of magnitude of a few hundred Kelvin, proba- 
bly somewhat overestimates the spin gap since it ignores 
correlations and disorder effects. 

The chiral spin-polarized edge states offer a unique set- 
ting to study spin transport. In particular, the spin-split 
state v = may be used to generate and detect spin- 
polarized currents. This spin transport regime seems at- 
tractive due to the large bulk gap and high stability of 



B 



FIG. 5: A Hall bar at v — can be used to generate and 
detect spin currents. Blue and red lines represent edge cur- 
rents with up and down spins. Contacts 1 and 4 are source 
and drain, which may be used to inject spin polarized cur- 
rent. Contacts 2, 3 are voltage probes with full spin mixing. 
The measured Hall voltage is directly related to spin current 
flowing in the system. An asymmetry between the upper and 
lower edges, e.g., introduced by removing voltage probe 3 or 
by gating, creates spin filtering effect: an unpolarized current 
injected from source 1 induces a spin-polarized current flow- 
ing into drain 4. Hall probes 5 and 6 downstream can serve 
as detectors of spin currents. 



the edge states. Moreover, increased quality of samples 
should allow existence of spin polarized edge states even 
at relatively low magnetic fields. 

Since the purpose of this section is mostly illustrative, 
we will keep our discussion as simple as possible. In 
particular, we shall ignore transport in the bulk, leaving 
the discussion of its role for Sec. IV. We also first neglect 
spin flip backscattering between edge states within one 
edge. Estimates of the spin flip rate wil be given below, 
Eq. lfMj) . A general approach, based on the Landauer- 
Biittiker formalism |24| , which can be used to calculate 
spin and charge currents at the edge for any configuration 
of current and voltage leads, was presented in Ref. 0. 

In this approach, transport is described by a scattering 
matrix [24J , with the edge states playing the role of scat- 
tering channels, and the reservoirs supplying in-states 
and absorbing out-states. Current in each mode is de- 

2 

scribed by the relation I^rn — Tr^tCi)' wnere VT(l-) * s ^ ne 
reservoir chemical potential for given spin projection. We 
consider the Hall bar geometry with four contacts 1-4 (see 
Fig©, where the contacts 1 and 4 serve as current source 
and drain. For these two contacts we do not assume spin 
mixing, so that the injected and drained current may 
be spin polarized. The contacts 2, 3 are voltage probes, 
which means that they do not drain current from the 
system. Furthermore, we assume that the probes pro- 
vide full spin mixing, i.e. chemical potentials of outgoing 
spin-up and spin-down electrons are equal. 

The simplest situation arises when unpolarized current 
is injected through contact 1. Then the up- and down- 
spins spatially separate in a symmetric way, flowing along 
the opposite edges of the bar. This can be interpreted 
as circulating spin current, and described as spin-Hall 
effect with quantized spin conductance Uxy"^ — e 2 /h. 
No electric voltage will be induced between the voltage 
probes 2, 3 in this case (zero charge-Hall effect). 
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This device can be used as a detector of spin polarized 
current, made possible by the reciprocal of the spin Hall 
effect, in which the electric Hall voltage is directly related 
to spin rather than charge current. Suppose the up-spin 
and down-spin electrons, injected through contact 1, have 
unequal chemical potentials, ip-\ ^ tp^. Then the cur- 

2 

rents flowing into the probes 2 and 3, = x^TCi)' 

after equilibration and spin mixing in the probes, in- 
duce voltages V2(3) = <^q)/2. The resulting Hall volt- 
age V xy = (ipi — ¥>t)/2 is directly proportional to spin 
current. At the same time, an unpolarized current (for 
which ip i = ipi) flows symmetrically in the upper and 
lower edges without generating Hall voltage. 

Spin transport at v = also allows to realize spin filter. 
Suppose that the upper and lower edges of the device in 
Fig|5] are made asymmetric, which can be achieved, for 
example, simply by removing probe 3. Then wc inject 
unpolarized current 21 through contact 1. The injected 
current will be distributed equally between the upper and 
lower edges in cross section A. In cross section B, how- 
ever, the net current will be spin polarized due to spin 
mixing in probe 2. The down-spin current reaching the 
drain in the upper edge equals 1/2 while the up-spin cur- 
rent in the lower edge is i\ Therefore, the total drained 
current becomes spin polarized. The spin polarized cur- 
rent can be fed into another system (see Figj5j) , where it 
can be detected using Hall probes 5 and 6 as discussed 
above. 

More complicated circuits can be assembled which gen- 
erate spin currents and detect them elsewhere. Note that 
the important principle is that as long as backscattering 
is not allowed, the edge current can travel long distances 
and the circuit is nonlocal, just as in the integer QHE 24]. 
In this case the current-voltage relationship is obtained 
by solving the circuit equations as described in Ref.[9(. 
The spatial scale of nonlocality is controlled by spin re- 
laxation which can be due to spin-orbit interaction or due 
to magnetic impurities near graphene edge. 

For simplicity, here we limit the discussion to the ef- 
fects of spin-orbit. There are two main spin-orbit terms 
in the graphene Hamiltonian [25| , the so-called intrin- 
sic and Rashba interaction, given by 



Hso = ^soVzT z s z , H R = A R (a x T z s y - a y s x ) , (23) 

where Pauli matrices <7i act in the sublattice space (Dirac 
spinor), while act in the valley space, and Sj represent 
physical spin. Estimates from band calculations [25[ give 
Ar « 0.1 K and a negligibly small Aso ~ 6mK. 

To estimate the backscattering rate due to the spin- 
orbit interaction (|23j) . we note that for an ideal atom- 
ically sharp edge the spin-orbit would couple the left 
and right states with the same momentum, opening a 
minigap at branch crossing: e = ±\/ e 2 + Ar 2 . How- 
ever, this momentum-conserving interaction alone can- 
not backscatter edge states, and we need to take disorder 



into account. Edges of graphite monolayers have been 
imaged using STM probes [13|, |l4|, where it was found 
that typically edge disorder can be viewed as patches of 
missing atoms of characteristic size d ~ 1 nm. Taking 
into account the left-right branch mixing by spin-orbit 
i?R, characterized by small mixing ratio of Ar/|£| <C 1 
away from branch crossing, we obtain an estimate of the 
backscattering mean free path: 



t{e) ~ {e/X R ) 2 (e B /d) 2 d, 



> 



Ar, 



(24) 



which gives I ~ 10 /zm for typical e ~ 10 K. The factor 
(£s/d) 2 > 1 accounts for the magnetic field dependence 
of disorder matrix elements. 

The quadratic energy dependence in (|24[) . with spin 
flip rate having a sharp peak near branch crossing, sug- 
gests @ the possibility to control backscattering using 
local gate. By tuning local chemical potential to and 
from the branch crossing, where the spin flip rate has a 
sharp peak, Eq. (j2~4"|) , we can induce or suppress backscat- 
tering in a controlled way. Spin filtering is achieved by 
controlling local gates on opposite sides of the Hall bar 
asymmetrically. 



IV. EDGE AND BULK TRANSPORT AT v = 0. 

Spin flip backscattering (|24|) can be incorporated in the 
edge transport model, described by coupled equations for 
particle density in the two spin-polarized modes: 



d t ni + d x ip x = 7(^2 - Wi) 
d t n 2 - d x (p 2 = 7(^1 - 92), 



(25) 



where 7™ 1 = £ is the backscattering mean free path 
taken for e at the Fermi level, and v± t 2 are compressibil- 
ities of the modes. (For brevity, we use 1 and 2 instead 
of | and J..) In writing Eqs. (l25p we implicitly assume 
that fast energy relaxation maintains local equilibrium 
of each of the modes, which is consistent with metallic 
temperature dependence of transport coefficients [17| • 

In a stationary state, Eqs. ([25|) have an integral I = 
fi — <^2 which expresses current conservation at the edge. 
[In this section we use the units of e 2 /h = 1 .] The general 
solution in the stationary current-carrying state is 



Vi,2 (as) = ft* - £ = ~fl 



(26) 



Taking into account that / is the current in one edge, we 
calculate the total current as 



I = 21 = -£ 

1 



(27) 



To describe the longitudinal resistance in the four- 
terminsl geometry, one must add potential drop on volt- 
age probes which gives R = ^(jL +1), where L is 
the distance between the probes. Comparing to the data 
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for p xx at v — we estimate till £ ~ 0.5 fim. This mean 
free path value, which is relatively small on the scale pre- 
dicted by Eq. (j2"4"|) , can be explained if spin flip processes 
are dominated by nonintrinsic effects, such as magnetic 
impurities localized near the edge. 

It is crucial that the edge transport model (|25|) treats 
both edges of a Hall bar in an identical way, thus pre- 
dicting zero Hall effect. In order to understand the 
observed density dependence of Hall coefficient 

us ma, 

which changes sign smoothly at v = without exhibit- 
ing a plateau, and of p xx which has a sharp peak at 
v = 0, we need to incorporate transport in the bulk in our 
model. In the full edge+bulk model, the density depen- 
dence of transport coefficients arises from bulk currents 
short-circuiting edge currents away from v — 0. This ex- 
plains, as we shall now see, the Hall effect, the peak of 
p xx , the resistance fluctuations near v = 0, as well as the 
behavior of a xx and a xy . 

We describe the transport problem in the bulk by the 
current-field relation, separately for each spin projection: 



-CTiVlpi, Oi 



& xy 
Uxy v xx 



1,2, 



where ip\.2 are electrochemical potentials for two spin 
states. We assume that the bulk conductivities a xx , as 
a function of density v, are peaked at the spin-split Lan- 
dau levels. For simplicity, here we ignore possible valley 
splitting, in which case the spin up and down Landau 
levels occur at v = ±1 around the Dirac point. As a 
simplest model, below we use Gaussians 



r(2) 



-AO+1) 2 



(28) 



with the parameter A describing the width of the levels. 

(X 2) 

The Hall conductivities a X y exhibit plateaus on either 



side of the peak in a-hx ■ The dependence of a' X y^ ) on v 
can be modeled with the help of the semicircle relation 
<?xy 2 ^ (axy T 2) + (cii' 2 ^) 2 = which often provides a 
good description of conventional QHE systems [26[ . 

The condition of charge continuity, Vji = 0, gives 
a 2d Laplace's equation for the potentials, V 2 -0i = 0. 
This equation must be solved together with the bound- 
ary conditions phenomenologically describing bulk-edge 
coupling: 



(1,2) 



n.j 2 



(29) 



where n is a normal vector to the boundary, and g{ip% — 
ipi) represents the edge-bulk leakage current density. 

Although a general solution of this problem can be 
given with the help of Fourier method, here we consider 
only the case when the potentials ipi (x) are varying slowly 
on the scale of the bar width w, which will suffice for our 
analysis of a homogeneous current flow. In this case, 
linearizing ipi(x) in the direction transverse to the bar, 



we can write Eqs. (l29]) for both edges of the bar as 

-a- xv d x ipi + a xx {ip v - ipi)/w = g(ipi - <pi) 
<y xy d x ipi> + a xx (ipi - ipi>)/w = g(ipi> - ipi>), 



(30) 



i = 1,2, where the primed and unprimed quantities de- 
note variables at opposite edges of the bar. 

Equations for the edge variables (fi are obtained by 
adding the bulk-edge leakage term to Eqs. (|2"5")) . giving 



(31) 



d x ipi = 7(^2 - <pi) + g(i>i - fi), 

-d x (p 2 = j(<pi - fa) + 3(^2 - ^2), 



along with a similar pair of equations for ipy , ip2> at the 
opposite edge. 

The solution of these eight equations, describing uni- 
form current, is of the form ipt = (p* — £x, ipi = ijj* — £x, 
etc., with the same linear part —Ex for all quantities. 
Using the algebraic structure of this linear system and 
the symmetry between the edges, we reduce the number 
of equations from eight to two. First, it is convenient 
to express the parameters 99* through ip* using Eqs. tpTLj) . 
which gives 



f2 



7+g ,/.* 

2-y+gVl 
27+9 % 



27+9 

7 / * 

27+^1 



r 2 + 



!!! 

27+9 



(32) 



Writing similar equations for the variables at the opposite 
edge to express ipi', <p2' through ipy, ip2', and substitut- 
ing the result in Eqs. (|30p . we obtain four equations for 
ipi and ipi' which have the form 



-d ( xy ] w£ = a ( xx \r v ~ i>*) + A(V>! " ^1) 

&xy w£ = (TxJ(lp* —1p*>) + \{lp2' ~ i'l') 



-a x 2 Jw£ = a x 1>(iP* 2 , - iP* 2 ) + \{iP*i ~ r 2 ) 
a x 2 y ] w£ = ai 2 x \iP* 2 - ip{) + \{iP* v - iP*,) 



r (2) 



(33) 



where the coefficients in this linear system are defined as 



£(1,2) 



,(1,2) 



dz 
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2 7 + g> 



27 + £ 



(34) 



~ (1 2) 

The quantities a X y represent the sum of the bulk and 
edge contributions to Hall conductivity for each spin. 

Symmetry between the edges allows to further reduce 
the number of independent variables. For that we add the 
first two equations to obtain ip\ + ip*, = ip 2 + 1P2' ■ Also 
we note that all potentials can be changed by the same 
constant that can be chosen so that the new quantities 
ip* and ipi satisfy ip* = -ip*,. After that Eqs.d25]) yield 



~(1) c n (1) /* 

aiywc = Zaixtpl 
a x 2 yw£ = 2aiiVl 



A(^ 2 * 



1$) 



(35) 



These two equations can be solved to find ip* 2 ■ 

Now we can find the current as a sum of the edge and 
bulk contributions, / = / c d ec + ^1 



'bulk 



where 



<pi — tp2 + <P2' ~ <Pv = ^{Vi - V2) 
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FIG. 6: Density dependence of transport coefficients 
pxx = Jw/2, p xy= ^/ 2 and = p xy /(p xy 2 + p X x 2 ), 

G X y = Pxy/(pxy 2 + pxx 2 ), obtained from the edge transport 
model (|31[) augmented with bulk conductivity, Eqs. (|30[l (see 
Eqs, (|36|l . (|37|l and text). Parameter values: A = 6, "fw — 5. 
Note the peak in p xx , the smooth behavior of p xy near v = 0, 
a quasi-plateau in G xy , and a double-peak structure in G„. 



can be discerned from the double peak structure in G xx 
in Fig(6j The peaks correspond to the bulk Landau level 
contributions, Eq. (f!?5)) . whereas the part of G xx between 
the peaks, exceeding the superposition of two Gaussians, 
Eq. (f!?5|) , is the edge contribution. The Hall resistance 
p xy is nonzero due to imbalance in a xy for opposite 
spin polarizations away from v — 0. Interestingly, p xy in 
Fig|6] exhibits no plateau, while G xy calculated from p xy 
and p xx displays an under-developed plateau-like feature. 
Overall, this behavior resembles that of the experimen- 
tally measured transport coefficients [IF 



17] 



Another notable feature of the measured p xy and p xx 
is enhanced fluctuations near zero v. These fluctuations 
are found to be strong in Ref . (l6| . where p xy changes 
sign several times near v = 0. They are also present, 
although are not as dramatic, in Ref. [171]. In the latter 
case, both p xy and p xx exhibit noisy behavior in the in- 
terval near v — comparable to the p xx peak width. As 
Ref. [13] points out, this behavior is consistent with the 
edge transport model. In the absence of bulk transport, 
the distribution of potential along the edge depends on 
the local backscattering rate 7(2), whereby Eq. (j2"fj|) is 
replaced by 



and 



/bulk = 4V (Vi - M + 42 + 4? ( v>2 - ) + 42 we 



After expressing <pt through ip, with the help of Eqs. 
and using the solution of Eqs. ([55]) . we obtain a relation 
I = 2£/7, where 



2 

1 



2 1 + g 



o (1) 

l J XX 



(2) 
Pxx 



. , (1) .(2) , (2 

AW {(T xy /<J XX - <Txv/<Ti: 

2 + \/aj& + A/4? 



.(36) 

The quantities plz' 2 '' are defined as p xx = 4x/(4y 2 + 

4 l 2 2 ). The quantity 7, Eq.([5B]). replaces 7 in Eq.([27p. In 

(1 2) 

the absence of bulk conductivity, a xx — > 0, we recover 
the result for pure edge transport, 7 = 7. 

The Hall voltage can be calculated from this solution 
as Vh = \{<£i + <-P2 — fv — f2'), where <fi, <fi> are variables 
at opposite edges. We obtain Vh = ££, where 



^ = 2w 



O'xy ^A + <T xx ^j + <j xy ^A + (TxxJ 



2a (1 V 2) 

XT u XX 



a4J 



a4J 



(37) 



This quantity vanishes at v — 0, since a. 



(i) 



a (1) 

u XX 



.(2) 



-a (2) 



and 



a xx at this point due to particle-hole symmetry. 
Transport coefficients, obtained from this model for 
typical parameter values, are displayed in Fig[6]which re- 
produces many of the key features of the data (see Fig.l 
in Ref. fl7l] ) . In particular, the peak in p xx is due to edge 
transport near v = 0. The suppression of p xx at finite v 
is due to the bulk conductivity short-circuiting the edge 
transport. The bulk and edge contributions to transport 



V?i,2(a0 = Vi,2(0) 



I / j(x')dx'. 
la 



Fluctuations of 7 arise due to its sensitivity to the local 
value of Fermi energy in the spin-orbit scattering model, 
Eq. (|24[) , and, similarly, for the magnetic impurity scat- 
tering mechanism. Assuming that the random part of 7 
is of a white noise character, we obtain strong fluctua- 
tions 8<p\fl{x) along the edge of magnitude that scales 
as a square root of the edge length. These fluctuations 
will contribute equally to the longitudinal and transverse 
voltage, since they are uncorrelated on the opposite sides 
of the Hall bar. The absence of fluctuations away from 
v = can be understood as a result of bulk conductivity 
short-circuiting the edge current, which will equilibrate 
potentials on the opposite sides of the Hall bar. 

The above discussion summarizes the results drawn 
from an attempt to model quantum Hall transport in 
graphene at v = by counter-circulating edge states. By 
taking into account backscattering within one edge as 
well as conduction in the bulk which short-circuits edge 
transport away from the neutrality point, this model ac- 
counts for the observed behavior of transport coefficients. 
Still, since no direct evidence for spin polarization has yet 
been found, more experimental and theoretical work will 
be needed to confirm the chiral spin-polarized edge pic- 
ture of the v — state. If proven to exist in graphene, 
these states will provide a unique setting to study spin 
transport as well as other interesting phenomena. 
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